Adiabatic and isothermal sound waves: the case of supercritical nitrogen 
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The acoustic sound dispersion of nitrogen in its liquid and supercritical phases has been studied 
by Inelastic X-Ray Scattering. Approaching supercritical conditions, the gradual disappearance of 
the positive sound dispersion, characteristic of the low temperature liquid, is observed. In the super- 
critical state, evidence for a crossover between adiabatic and isothermal sound propagation regimes 
is inferred by an analysis of the dynamic structure factor based on generalized hydrodynamics. 
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The development of Inelastic X-ray Scattering (IX S) 
EJ LH LI has opened up new opportunities in 
the study of the high-frequency dynamics of liquids at 
nanometer distances and picosecond timescales. Various 
IXS investigations unveiled that, in the proper thermo- 
dynamical conditions, structural relaxation processes oc- 
cur in the THz frequency range. These processes mani- 
fest themselves by an upward bending of the sound dis- 
persion relation, which deviate from the low-frequency 
linear behavior dictated by the adiabatic sound veloc- 
ity. Accordingly, the sound velocity increases from the 
adiabatic value, c*, to the "infinite" freq uency value, 
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If 0_l(Q) is the frequency of the probed acoustic exci- 
tation with wave vector Q, and t q is the characteris- 
tic relaxation time of the structural relaxation process, 
the system exhibits a liquid-like (viscous) behavior if the 
condition Ql,(Q)tc <C 1 is fulfilled: in this limit the dis- 
persion relation is Ql(Q) — c s Q. In the opposite limit, 
the elastic one, where Vli,(Q)T a S> 1, structural rear- 
rangements are too slow to efficiently dissipate the en- 
ergy of the acoustic waves: in this limit Q,l(Q) = c^Q. 
A schematic picture of this visco— elastic variation of the 
sound velocity form c s — * Coo is shown in Fig. la (full 
thick line), where the cross-over condition VLl(Q) = 1/t q 
is indicated by the vertical arrow, and the adiabatic and 
infinite linear dispersion relations are indicated by the 
the full and dashed lines respectively. Here the dash-dot 
line represents the value of 1/t q , which is here assumed 
to be Q-independcnt. 

The structural relaxation process, in principle, is not 
the only relaxation process that could affect the high 
frequency dynamics of a liquid. On general grounds, 
also other physical processes, as for example thermal 
diffusion, can modify the behavior of sound waves. In 
this case, given a thermal diffusion coefficient Dt and 
a thermal relaxation time tt — 1/DtQ 2 , a modifica- 
tion of the sound velocity takes place around a cross- 
over frequency given by Ql(Q) ~ DtQ 2 - In this case, 
for DtQ 2 <C ^l(Q) the thermal diffusion is much 
slower than the period of acoustic waves, which there- 
fore propagate adiabatically, i.e. without any thermal 
exchange with the local environment. In the opposite 
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FIG. 1: Schematic representation of the various sound regimes 
highlighting the two cross-over conditions. Upper panel: vis- 
coelastic transition taking place at Q.l — 1/t; middle panel: 
isothermal transition taking place at Ql = DtQ 2 ; bottom 
panel: realistic non linear dispersions under favorable exper- 
imental conditions for Q,l <C 1/t and Q.l ~ DtQ 2 ■ Verti- 
cal arrows indicate the crossover regions as discussed in the 
text. The isothermal (dotted line), adiabatic (solid line), in- 
finite frequency (dashed line) and apparent (thick solid line) 
sound dispersions are reported together with the relaxation fre- 
quency 1/t (dashed- dotted line) and thermal frequency DtQ 2 
(dashed- double dotted line). 



limit, DtQ 2 3> £Il(Q), the thermalization of the acous- 
tic wave is instantaneous, and the effective propaga- 
tion mechanism becomes isothermal with sound velocity 
ct = ^~ x l 2 c s (7 is the specific heat ratio). One there- 
fore expects to observe a transition between the adiabatic 
and the isothermal regimes when £Il(Q) = DtQ 2 - This 
case is schematically presented in Figure lb, where lin- 
ear lines represent the adiabatic (full line) and isothermal 
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(dotted line) sound dispersion, and the dash-dot-dot line 
indicates the Q dependence of 1/t t (Q). 

Contrary to the visco-elastic transition, which has been 
experimentally documented in many liquids 0, H, IS IS 
IS 13 IH El ES 13 13 ES El EE no firm evidence of 
the adiabatic to isothermal transition is available. This is 
mainly due to the following reasons: i) The amplitude of 
the dispersive effect between c s and c-r is proportional to 
7 — 1 which, in most liquids, is close to zero, ii) The effect 
is often completely masked by the competing visco-elastic 
transition, iii) The crossover occurs typically at large Q, 
often above the position of the first sharp diffraction peak 
(FSDP) of the static structure factor, S(Q). Here, how- 
ever, the finite-Q generalization of the isothermal sound 
speed ct{Q) - which is proportional to 1/S(Q) 1 ^ 2 - [l9|. 
becomes too small to be reliably extracted from the mea- 
sured S(Q,u>) spectra. 

These obstacles to the observation of the adiabatic 
to isothermal sound propagation transition can be re- 
duced close to supercritical conditions for the following 
reasons: i) the competing positive dispersion induced 
by the viscoelastic transition is expected to be much 
weaker, ii) The crossover condition for the viscoelastic 
transition is expected to move at higher Q, i.e. higher 
frequency, since the structural rearrangements occur on 
shorter time scales, iii) Moreover, the intensity of the 
FSDP is much weaker, and consequently ct(Q) does not 
display any longer a pronounced minimum in the vicinity 
of the FSDP. iv) Finally, the sound speed decreases on 
approaching the critical point, and therefore the transi- 
tion will move towards smaller Q values. This scenario is 
schematically pictured in Figure lc, where the combined 
effects of the structural and thermal relaxation processes 
are model within the framework of the molecular hydro- 
dynamics The values of r Q and tt(Q) are chosen 
such that the transition is actually determined by the 
thermal relaxation process, as indicated by the arrow. 

This letter, guided by the above considerations, re- 
ports an IXS study of liquid and supercritical nitrogen. 
Nitrogen has been chosen because the two thermody- 
namic regimes (liquid and supercritical) can be easily ob- 
tained within an experimentally accessible pressure and 
temperature range, and because of its favorable inelas- 
tic scattering cross section. Figure 2 shows the portion 
of the explored thermodynamic plane. The inset reports 
the expected Q- values, Q*(T), for which the crossover 
from the adiabatic to the isothermal regime should oc- 
cur: Q*(T) — c s /Dt- The values of c s and Dt were 
obtained from the nitrogen Equation-of-State (EoS) [2fT |. 

The experiment was performed on ID28 at the Euro- 
pean Synchrotron Radiation Facility (ESRF). The beam- 
line was set up with an instrumental resolution function 
of l.5meV full- width-half-maximum (FWHM). The sam- 
ple was embedded in a large volume, high pressure cell 
kept in thermal contact with the cold finger of a cryo- 
stat. The optical windows were two 1 mm thick dia- 
monds disks, and the sample length amounted to f mm. 
The pressure stability was better than 5 bar over the ac- 
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FIG. 2: Phase diagram of nitrogen. The investigated ther- 
modynamic points are reported as full circles. The solid line 
represents the coexistence curve, ending at the critical point 
CP (T c =126.2 K, P c = 34 bar, p c = 0.313 g/cm3). In the in- 
set we report the expected Q-values for the isothermal transi- 
tion for the sound propagation, as estimated from the nitrogen 
equation- of- state using Q — c s /Dt- 



quisition time of a typical spectrum. IXS spectra were 
recorded following an isobaric path at 400 bar from sub- 
critical to supercritical conditions (full dots in figure 2). 
Cell contribution as determined by empty cell measure- 
ment yielded a negligible contribution to the signal for 
all exploited scattering geometries. Multiple scattering 
has been estimated to be negligible. The recorded signal 
is therefore proportional to the dynamic structure fac- 
tor, S(Q, uj), convoluted with the instrumental resolution 
function. 

The data analysis has been performed using a line- 
shape model derived within the framework of the mem- 
ory function formalism [Tsj . In this model the dynamic 
structure factor is written as: 



KbT 



[n(uj) + 1] Im[u 2 
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q c T 



iujm Q (uj)] 1 
(1) 

where n(oj) is the Bose factor and ttiq^uj) is the Fourier 
transform of the time dependent memory function 
Il9j . which has been approximated by: 



m Q (t) = 4(Q)[7(Q) ~ l]Q 2 e- DT(Q)Q 4 + 

+[cUQ) - cl(Q)]Q 2 ^ t/TaiQ) + A 2 (W) 



(2) 



Following the prescriptions of molecular hydrodynam- 
ics, all the thermodynamic quantities entering Eq. [3 are 
assumed as Q dependent extensions of their respective 
macroscopic counterparts. 

The 8 function in Eq[21 accounts for the fast dynam- 
ics which induces an instantaneous time decay of vtlq (t) 
H |S] • The second exponential term of Eq|21 accounts for 
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FIG. 3: Selection of IXS spectra of liquid nitrogen at P = 400 
bar and two temperatures (87 K and 190 K) at momentum 
transfers as indicated in the individual panels. The experi- 
mental data (full circles) are reported together with their best 
fit lineshapes (solid lines). 



the visco — elastic transition that disperses the sound 
velocity, with increasing Q, from its low to its high fre- 
quency limiting value, c s (Q) and Coo(Q) respectively. Fi- 
nally, the first term in Eq. [2] accounts for the thermal 
diffusion process. In the low-Q limit the thermal part of 
the memory function yields a central peak in the spec- 
tral density, whose integrated intensity is proportional 
to j(Q) — 1, and its half-width is given by Dt{Q)Q 2 ■ 
Consequently the proposed model predicts the "classi- 
cal" hydrodynamics description at low Q. 

The experimental data are described by the convolu- 
tion of the model reported in Eq.sQ] an d[21with the ex- 
perimentally determined instrument resolution function. 
In the fitting procedure, the parameters r a (Q), Coo(Q), 
A 2 (Q) and an overall intensity factor A are left as free 
parameters. Owing to the lack of experimental or com- 
putational results on the Q-dependence of 7 and Dt, 
these variables have been fixed to their thermodynamic 
{Q — > 0) values, as derived from the EoS of nitrogen [2oT |. 
The parameters ct{Q) and cJQ) have been fixed to the 
values calculated through |19|: 



ct(Q)=1- 1/2 c s (Q) = { 



K B T 
MS(QY 



(3) 



where Kb and M are the Boltzmann constant and 
the molecular mass of nitrogen, respectively. The static 
structure factor S(Q) was determined experimentally by 
recording the energy-integrated scatte ring , corrected for 
the known nitrogen form factor, f(Q) [21|, as well as for 
all geometrical artifacts. In order to put S(Q) on an 



absolute scale, we have normalized the measured inten- 
sity to fulfill the compressibility constraint in the Q = 
limit. Finally, the apparent sound velocity of acoustic 
excitations, cl{Q), has been determined using the rela- 
tion Cl(Q) = Ql(Q)/Q, where Qt(Q) is the frequency 
value of the maximum of the longitudinal current spec- 
tra, uj 2 S(Q,uj)/Q 2 . 

Selected IXS spectra and their best fits are shown in 
Fig. 3. Each spectrum typically covers an energy transfer 
range of ±40 meV and was collected in the 2 to 14 nmT 1 
Q range. The logarithmic plot emphasizes the overall 
good agreement between experimental and model line 
shapes, even in the tails of the spectra. 

The comparisons between the apparent dispersion re- 
lation, Ql{Q) vs. Q (open circles), and both isothermal 
(dotted line) and adiabatic (full line) ones, are shown in 
Fig. 4 at the four different temperatures 87, 128, 171, 
and 191 K. In the same figure we also report the infinite 
frequency dispersion relation (full dots) and the inverse of 
the structural relaxation time l/r Q (Q) (dashed-dot line) 
as derived from the best fit result. The inverse of the 
thermal relaxation time, DtQ 2 , is also reported as de- 
rived from the nitrogen EoS (dash-dot-dot line). The re- 
sults of Fig. 4a, corresponding to the lowest investigated 
temperature (87 K), show that the apparent dispersion is 
systematically higher than the adiabatic one. Although 
the Coo -limit is not completely reached, in this thermo- 
dynamic point the high frequency dynamics of the liq- 
uid is mostly influenced by the structural relaxation pro- 
cess. This influence is considerably weaker at T = 128 K 
(Panel b) due to the higher value of l/r a (Q) which pre- 
vents even-more to reach the viscoelastic crossover con- 
ditions. On further increase of the temperature (Panel 
c and d), l/r a (Q) moves definitely out of the probed 
frequency range and any hint of positive dispersion on 
the sound velocity disappears. In fact, with increasing 
temperature the apparent dispersion relation initially ap- 
proaches the adiabatic sound dispersion, but at the high- 
est temperature is observed to go below the adiabatic 
dispersion and to approach the isothermal one: this is 
the onset of the expected transition between adiabatic 
and isothermal regimes. This effect is most clearly seen 
for the three highest Q points between 9.5 and 14 nmT 1 . 

The arrows in Figure 4 (panels c and d) indicate the Q 
points where one would expect the adiabatic-isothermal 
transition (Q,t,(Q) = DtQ 2 , left arrows) and where the 
transition actually occurs (right arrows): a Q value 
definitively higher than the one estimated by the left 
arrow. This difference between expected and observed 
transition Q points is an indication of the Q-dependence 
of the parameter Dt, which we neglected so far. This be- 
havior is in agreement with the expectations suggeste d by 
numerical simulations of various molecular liquids |22Ll23j 
pointing to a monotonic decrease of Dt(Q) with increas- 
ing Q. As a consequence, the position of the crossover 
is shifted towards higher momentum transfers. Even 
though the transition to the isothermal regime is not fully 
completed, its occurrence is inferred from this system- 
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atic trend. Furthermore, it is worthwhile recalling that 
the present results have been obtained without impos- 
ing any constraint on the line shape fitting parameters. 
Moreover, the values of the isothermal sound dispersions 
have been determined independently through the mea- 
sured S(Q). 

In conclusion our observations can be summarized as 
follows: i) the apparent sound branch of nitrogen show 
clear signatures of a positive sound dispersion in the liq- 
uid phase, thus witnessing the presence of a structural re- 
laxation process, ii) The positive sound dispersion grad- 
ually disappears while reaching supercritical conditions, 
most likely due to the shift of l/r a {Q) above the Ql(Q) 
values, iii) At the highest investigated temperatures, the 
transition from the adiabatic to the isothermal sound 
propagation is observed. 

These findings indicate that, in the liquid phase, the 
high frequency dynamic of a simple, non-associated fluid, 
such as nitrogen, is mainly ruled by the structural relax- 
ation, as manifested by the positive sound dispersion. In 
contrast, in the supercritical phase, the role of the struc- 
tural relaxation becomes negligible and the most evident 
dispersive effect is the bending down of the apparent dis- 
persion, from its adiabatic value to the isothermal one. 
In this regime, the thermal relaxation process rules the 
high frequency dynamics. 



FIG. 4: Dispersion curves of the longitudinal acoustic mode 
of nitrogen at 87, 128, 171 and 190 K. The full and dotted 
lines represent the adiabatic and isothermal dispersions, de- 
rived from S(Q) measurements, as discussed in the text; the 
open circles indicate the maxima of longitudinal current spec- 
tra. The dash-dotted and dash-double dotted lines represent 
1/t(Q) and DtQ 2 , respectively; the former has been evaluated 
interpolating the results of the performed lineshape analysis, 
which allowed also to extract the values of infinite-frequency 
sound dispersion, reported as full circles. 
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